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Abstract. We explore a nonlocal connection between certain linear and nonlinear 
ordinary differential equations (ODEs), representing physically important oscillator 
systems, and identify a class of integrable nonlinear ODEs of any order. We also devise 
a method to derive explicit general solutions of the nonlinear ODEs. Interestingly, 
many well known integrable models can be accommodated into our scheme and our 
procedure thereby provides further understanding of these models. 



1. Introduction 

In a recent paper, we have shown that the modified Emden type equation with additional 
linear forcing, 

x + 3kxx + k 2 x 3 + \x = 0, (1) 

where over dot denotes differentiation with respect to t and k and A are arbitrary 
parameters, exhibits certain unusual nonlinear dynamical properties [1]. For a particular 
sign of the control parameter, namely, A > 0, the frequency of oscillations of the 
nonlinear oscillator Q is completely independent of amplitude and remains the same as 
that of the linear harmonic oscillator, thereby showing that the amplitude dependence of 
frequency is not necessarily a fundamental property of nonlinear dynamical phenomena 
[1]. In this case, namely, A > 0, the system admits the explicit sinusoidal periodic 
solution 

A sin (tut + 5) lo /— 

x(t) = -j— ± r—t- -, 0<A<-, u = VX, 2 

where A and 6 are arbitrary constants. The system exhibits certain other unusual 
properties also. For more details on the dynamics of this equation one may refer to [1]. 

In the same paper, we have also noted that the solution (J2J) can be derived 
unambiguously by introducing a nonlocal transformation, 

U = xe k Io (3) 
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and transforming equation (Q) to the linear harmonic oscillator equation, 

U + XU = Q. (4) 

For A > 0, the solution is U = Asm(ut + 5), where A and 5 are arbitrary constants and 
the frequency, u = y/X, is independent of the amplitude. Using (JBJ) one finds that 

Integrating Riccati type equation (jSJ) and absorbing the integration constant with the 
existing ones, one can obtain the solution (J2J). We mention here that the solution (J2J) 
can be derived in a number of ways, for more details one may again refer to [1] and [2]. 
One may also note that in the case A = the transformation (jSJ) modifies equation (0) to 
the free particle equation. The linearization and integrability properties of equation ([Tjl 
with A = has received a considerable attention in the nonlinear dynamics literature, 
see for example, Refs. [3-6]. 

Now the question that naturally arises is whether the nonlocal transformation 
which connects the linear and nonlinear oscillators, namely, Eqs. ((IJ and (J3J) is a rare 
one or there exists a wider class of nonlinear dynamical systems that are connected 
with linear oscillators in a hidden way? Our studies reveal that there exists a class of 
linear oscillators that are connected with nonlinear oscillator equations through nonlocal 
transformations. For example, let us consider a damped linear harmonic oscillator, 
instead of the undamped linear harmonic oscillator, that is, 

U + Cl U + c 2 U = 0, (6) 

where c\ and c 2 are arbitrary parameters, and consider a more general form of nonlocal 
transformation, that is, 

U = x(t)efi f WW i (7) 

where f(x(t)) is an arbitrary function of x. Now substituting (J2J) into (JHJ) we get a 
nonlinear ODE of the form 

x + (2f + xf + cj)x + (f 2 + af + c 2 )x = 0, (8) 

where prime denotes differentiation with respect to x. Interestingly, one can now see 
that for certain specific forms of the function /, the associated nonlinear ODEs become 
the well known and well studied models in the current nonlinear dynamics literature. 
To cite a few, first let us consider the case f(x) — kx in (|7jl. The associated nonlinear 
ODE, 

x + (c± + 3kx)x + k 2 x 3 + c\kx 2 + c 2 x = 0, (9) 

is nothing but the generalized modified Emden type equation (MEE) whose integrability 
properties have been studied in detail in Refs. [3-10]. It is one of the linearizable second 
order nonlinear ODEs which admits eight parameter Lie point symmetries [4]. 

Choosing f(x) = kx 2 in (jZJ), one gets the generalized force- free Duffing- van der Pol 
oscillator (DVP) equation 

x + (ct + Akx 2 )x + k 2 x 5 + kc!X 3 + c 2 x = 0, (10) 
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as its nonlinear counterpart. More details on the mathematical aspects and the 
underlying dynamics of this equation can be found in Refs. [9-13]. On the other hand 
fixing f(x) = Jj, one gets another integrable nonlinear second order ODE of the form [9] 

x + C\i, + c 2 x H 1 — - = 0. (11) 

The examples illustrated above clearly demonstrate that one can systematically classify 
a class of integrable nonlinear ODEs through this way. 

Motivated by these observations, in this paper, we carry out a detailed investigation 
on the nonlocal connection between certain linear ODEs and their nonlinear counterparts 
of any order. To make our study a systematic one we start our investigations at the 
level of second order and end with nth order equations. Besides constructing integrable 
nonlinear equations we also point out the importance of these nonlinear equations and 
obtain their solutions and thereby bring out the significance of these systems. 

The plan of the paper is as follows. In Sec. 2, we make a detailed study on the 
nonlocal connection between damped linear harmonic oscillator equation and certain 
nonlinear oscillator equations of second order and describe a method of deriving general 
solutions of the nonlinear ODEs. We extend this procedure to third and nth order 
ODEs in Sees. 3 and 4, respectively. We present a general theory which connects the 
variable coefficient linear ODEs with nonlinear ODEs and the method of finding the 
general solution in Sec. 5. We present our conclusion in Sec. 6. 



2. Second order ODEs 

To begin with let us consider a linear second order ODE of the form (jHJ), whose general 
solution is U = a(t), where a(t) is a known function. Now we consider a nonlocal 
transformation of the form 

u = x n e s t m , )* m +iV))dt> i (12) 

where n and m are constants and (3(t) and j(t) are arbitrary functions of t, and substitute 
it into (jnj) so that the latter becomes a nonlinear second order ODE of the form 

x+(n- I)— + ^x 2m+1 + h(t, x)x + b 2 (t)x m+1 + b 3 (t)x = 0, (13) 



x n 



where 



bi(t, x) = - (2nj + nci + (m + 2n)(3x m ^J . 
b 2 (t) = -($ + 2 1 f3 + c 1 (3 



h{t) = ^7 + 7 2 + 7Ci + c 2 ). (14) 

Obviously the transformation ()12j) corresponds to the special case, f(x) = 
/3(t)x m , j(t) = 0, in the transformation (J7J) but with the prefactor of the exponential 
function taken as x 11 , where n is arbitrary. We have been motivated to choose this form 
due to the examples illustrated in the introduction. In all these cases, one might have 
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noted that the function f(x) is a simple polynomial one. As a consequence, we would 
like to generalize this form and see the outcome. Such a generalization is indicated 
in Sec. 5. However, here we note that the form (|T3|) itself yields several important 
nonlinear ODEs as see below. 

Now the question is how to construct the general solution for the nonlinear equation 
(fT5Jl from the known solution of ©• This can be done by using the identity, 

U 



nx 



+ P(t)x m + 1 (t), 



(15) 



derived from equation (|12jl. Since U = a(t), where a(t) is a known function, the solution 
of the second order linear ODE equation (|15|) can be brought to the form 



x = (a(t) — 7(£) 



x 



-x 



m+1 



where a = 
the form 



(16) 

n n 

Solving equation (jTflj) we get the general solution for the equation (|T3j) in 



x (t) = e ^ Jo MVrC*'))*' \c +- f ( /3(t)e% Jo W')-7(*")K') dt 



n Jo 



(17) 



Note that the equation (JT7J) contains three arbitrary constants (two in a(t), since 
a(t) is a general solution of the second order equation (|B]L and another constant C). 
However, one can absorb one arbitrary constant with the other two and rewrite the 
solution in such a way that it contains only two arbitrary constants, as we see below. 

Let us assume that the solution of © be written in the form 

U(t) = a(t) = he rnit + I 2 e m2t , (18) 

where Ij's, i — 1,2, are integration constants and m^'s, i = 1,2, are the roots of the 
auxiliary equation associated with the ODE (JHj). From (JTSJl we get 

f/ _ mj^" 11 * + m 2 e m2t 

U ~ : 



a(t) 



a 
a 



J 1 grn, 1 t _|_ e m2t 



(19) 



where l\ = j^. Integrating both sides of the equation (fTH|) we get 



a{t')dt' = log{he mit + e 



log 



.a(t). 



(20) 



Substituting the expression ([2D)) into equation (|T7jl . we get a general solution for the 
equation ()13|) in the form 



x(t) 



' / 2 



X 



n Jo V i9 



)- e ^/o t '^") rf '"') ( it' 



(21) 



Now the general solution (|21|) contains only two arbitrary constants, namely, I\ and 
C. In the following, we consider certain important sub-cases of the linear oscillator 
and discuss their nonlinear counterparts and their physical significance. 
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2.1. Case (i) c\ = and c 2 = 

In this case the linear equation (JHJ) is nothing but the free particle equation. The 
associated nonlinear equation can be fixed from 1)13)1 by restricting c\ = and c 2 = 0. 
The general solution of the nonlinear equation turns out to be 



(22) 



C+- / [(5(t'){t' + h)~e 
n Jo V 

where ii and C are arbitrary constants. 

The interesting fact is that equation (|13)1 contains a family of important nonlinear 
ODEs. To mention one such example let us choose n — 1, j(t) = and /3(t) = k so 
that equation (fTH)l now becomes 

x+ (m + 2)fc£ m x + A;V m+1 = 0. (23) 

The invariance, integrability properties and direct linearization through generalized 
transformation of equation (j2l?)) can be found in Refs. [6,9, 10]. From our above results, 
the general solution of (|2*3*|) can be fixed easily from (J22j) in the form 

mk -(/i+rf, (24) 



x(t) = (t + h 



C + 



[m 



1 



which exactly coincides with the results obtained through other methods [6,9,10]. 



2.2. Case (ii) C\ = and c 2 = constant 

In this case, from equation ©, we have the linear harmonic oscillator equation at 
our hand. The nonlocal transformation (J12)l transforms the linear harmonic oscillator 
equation, U + c 2 U = 0, to the nonlinear ODE (fTB")) with c\ = 0. The general solution of 
dT3J) turns out to be 



x(t) 



cos(u;£ + 5) 



(25) 



x C + (/9(Oe- - ^ l{t " )dt " cosV(a;t' + S^dti] ™ , 

where u = ^fc~2 and 5 and C are arbitrary constants. 

Interestingly, one can deduce certain important nonlinear ODEs from (|13[) with 
Ci = 0. One such interesting case is n = 1, m — 1, 7(t) = and /?(t) = A;. In this case 
we get the equation (JTJ). The unusual dynamical properties exhibited by this equation 
has already been pointed out in the introduction. 

Now we choose n = 1, m = — 2, j(t) = and (3{t) = k in equation (JTSJ) so that the 
latter becomes 

k 2 

x + c 2 x + — = (26) 

x 6 

which is another important nonlinear ODE which arises in different areas of physics and 
has been studied in detail in Ref. [14,15]. The solution can be derived from ([25)1 in the 
form 



x(t) = cos(ut + 5) 



2k 

C tan(ut + S) 



C-2- 



(27) 
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Note that the solution defines a well defined harmonic periodic oscillation with period 
T = -2= which is the same as that of the unperturbed simple harmonic oscillator. 

This nonlinear oscillator again supports our argument that the amplitude dependence of 
frequency is not necessarily a fundamental property of nonlinear oscillations. 

On the other hand fixing n = 1, m — 2, j(t) = —k 2 and (3(t) = k\ in equation (jTBj) 
we get 

x + (Ahx 2 - 2k 2 )x + k\x b - 2k 1 k 2 x 3 + (c 2 + k\)x = 0, (28) 

which is nothing but the generalized force-free Duffmg-van der Pol nonlinear oscillator 
equation [9,11-13]. The general solution can again be fixed from (|2*K|) as 

x{t) = cos(u>t + 5) 



Ce 



-2k 2 t 



k x 



+ 2k 2 + 2(u 2 



k? 



-(&2 cos2(o>t + 5) + lo sin 2{ut + 5)) 



(29) 



"2 

which exactly coincides with the result obtained through other methods [9,11-13]. 

Finally, let us consider the case n = 1, m = q, j(t) = —k 2 and 0(t) = k\ in 
equation ()13jl . In this case we end up with the following nonlinear ODE 

x + {{q + 2)k lX q - 2k 2 )x + [{hx* - k 2 f + c 2 ]x = 0. (30) 

Equation ()30|) admits the following general solution 

i 



x(t) = cos(u;£ + 8)e 



k 2 t 



C + qh e qk2t ' cos q (ut' + 8) dt' 



(31) 



We mention here that for even values of q and positive values of k\ and k 2 equation (pffij) 
exhibits limit cycle oscillations [12]. 



2.3. Case (Hi) Ci, c 2 = constant 

In this case the nonlocal transformation connects the damped harmonic oscillator ODE 
with the nonlinear ODE (|13|) . From (|21j) . the general solution of equation (|13j) with c\ 
and c 2 taking constant values turns out to be 



x(t) 



2 



X 



n Jo V 



e (^)t' 



(32) 



where u> = y c\ — Ac 2 and I\ and C are arbitrary constants. 

In this case also one can identify several physically interesting nonlinear oscillator 
equations. For example, choosing n = 1, m = 1, j(t) = and (3{t) = k in equation 
()13|) we get an equation of the form (jUJ). The general solution can be written from (|32j) 
in the form 

2c 2 (h + e^) 



x(t) 



Ce { '^ L)t - k{c x - u))I x - k{c x + u)e ut 



(33) 



where to 



the Refs. [8,9]. 



4c 2 . Indeed our result exactly coincides with the result obtained in 
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The choice n — 1, m = 2, 7(i) = and f3(t) = k leads us to an equation of the 
form (fTU|) . The general solution of (fTU|) can be fixed easily from ()32j) in the form 

m = (-^ 2cMh+^? ^ (34) 

VC*e ( ^ _)t - kllci(ci -uj) - kci(ci + cj)e 2 ^* - %kl x c 2 e ui ' 
which indeed the general solution of the genaralized force-free DVP oscillator equation 
[9]. 

In this section, besides the general case, we considered three different linear 
equations, namely, free particle equation, linear harmonic oscillator and damped linear 
harmonic equations and demonstrated how they are related to the well known nonlinear 
models through nonlocal transformations. In the following section, we extend the theory 
to third order ODEs and study the outcome. 

3. Third order ODEs 

Let us consider a linear third order ODE of the form 

U +ci# + c 2 U + c 3 U = 0, (35) 

where C\, c 2 and C3 are arbitrary constants. The nonlocal transformation (j!2|) transforms 
(I35|) to the nonlinear ODE of the form 

■ ^ 

x +[3(n - 1)- + x)]x + (n - l)(n - 2)— + d 2 (t, x)x 2 

+d 3 (t, x)x + ^—x 3m+1 + d 4 (t)x 2m+1 + d 5 (t)x m+1 + d 6 (t)x = 0, (36) 

where 

di(t, x) = — (^3nj + nci + (3{m + 3n)x m ^j , 

d 2 (t,x) = -((m(m + 2n) + (n - l)(m + 3n)^J(3x m ~ l + n(n - + c 3 )x~ 
d 3 (t,x) = -^37^(m+2n) +/3(2ra + 3n) + c x (m + 2n)/3^z m 

+3/3 2 (m + n)x 2m + n(zj + 3-f 2 + 2ci7 + c 2 
d A {t) = ^3/3/3 + 37/? 2 + /5 2 ci) , 

d 5 (t) = i (p + 3(7/3 + 7/3) + 3/3 7 2 + (ci/3 + c 2 (3) + 2c l7 /3" 

d 6 (t) = - (7 + 377 + 7 3 + 01(7 + 7 2 ) + 7C 2 + c 3 ) . (37) 

To obtain the exact solution, x(t), of equation (}3l)j) we again apply the results 
of Eqs. (jl5j) - ([T7]) to the present case. Since we are using the same form of nonlocal 
transformation (jl2)l to transform the nonlinear ODE (j3(ij) to the linear ODE (|33|) . the 
general solution of (|36J) looks exactly the same form as ()17|). The only difference is that 
the general solution, in the present case, contains four integration constants (three of 
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them come out from the solution of the linear ODE (J35|) and the remaining one from 
solving (fT6*j) ). Now let us rewrite the general solution in such a way that it contains 
only three arbitrary constants by absorbing the fourth one. To do so let us follow the 
same procedure which we adopted in the case of second order ODEs. 
Let us write the general solution of equation in the form 

U(t) = a(t) = he mit + I 2 e m2t + he mz \ (38) 

where ij's, i = 1,2, 3, are integration constants and m^'s, i — 1, 2, 3, are the roots of the 
auxiliary equation associated with the differential equation (|33j) . From (}3"%)) we get 

. . a U mj x e m ^ + m 2 i 2 e m2t + m 3 e m3 * , , 

a(t) = - = — = , (39) 

a U ^e" 11 * + I 2 e m2t + e mst 

where Zj = ji, % — 1, 2. Integrating both sides of equation (}3"9"|) we get 

t a(t')dt' = \og(he mit + I 2 e m2t + e mst ) = log(^). (40) 
Jo h 

Substituting the expression (j4~Uj) into (fTTj). we get the general solution for the equation 
(J36|) in the form 

x(t) = (^e^i:^ dt '\c+- C (m(—^e^fo^'')dt']^. (41) 
^3 L n Jo \ J 3 / 

We note that the general solution (}4~Tj) now contains only three arbitrary constants, 

namely, Ii, I 2 and C. 

Since we have three arbitrary parameters in the linear ODE (|35j). one can consider 

certain specific sub-cases of physical interest. In the following, we discuss some of them. 

3.1. Case (i) a = 0,i = 1,2,3 

In this case, from equation ()35|) . we have a simple linear equation = which can 
be connected to the third order nonlinear ODE (J36|) (with ci, c 2 , c% = 0) through the 
nonlocal transformation (|12p. The general solution of ()36|) can be fixed easily from (|41j) 
in the form 

m f — 1 

„ \ n m P* / . 1 1 \ i . / / \ to 

(42) 



x 



C + f I i m (2 + ht ' + h ) " C " " ^ l{t " )dt ") dt ' 
where I\, I 2 and C are arbitrary constants. 

Equation ()3(j)l includes a class of important nonlinear systems. For example, 
choosing n — 1, m = 1, 7(£) = and (3{t) = k equation (|3^|) becomes 

x +4kxx + 3kx 2 + Qk 2 x 2 x + fcV = 0. (43) 

Equation (pf3"j) is a special case of the Chazy equation XII (with N = 2 and parametric 
restrictions A = 0, B = in Ref. [18]), which has been studied in detail in Refs. [16-22]. 
The general solution can be fixed easily from ()42j) in the form 

+ Lt + LIo 



kt 2 
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3.2. Case (ii) C3 = 0, C2 = and c\ = constant 

In this case, the linear equation (J35|) becomes [7 +cj.?7 = and the nonlocal 
transformation ()12j) transforms this ODE to the form (jMfijl with 02,03 = 0. In this 
case we have the general solution of the form 



x(t) 



x 



C 



m 
n Jo 



(3(t')(e- Clt ' +ht' + h) n e-^Io ^"W')dt' 

1, m 



(45) 

1, y(t) = and 



where and C are arbitrary constants. Choosing n 

Pit) = k equation (pjj) becomes 

x +(ci + Akx)x + 3kx 2 + 3k(a + 2kx)xx + (ci + kx)k 2 x 3 = 0. (46) 

The general solution of equation ()46|) can be fixed easily from ()45|) in the form 

e- Cl * + ht + I 2 



x(t) 



C + ^(2e-^ + c 1 I 1 t 2 + 2c 1 I 2 t) 



(47) 



3.3. Case (Hi) c\ = 0,C2 = and C3 = constant 

In this case, the linear ODE (|35|) assumes the form [/ +03^/ = and the nonlocal 
transformation (fT^j) transforms this equation, JJ +C3U = 0, to the nonlinear ODE 
with ci, C2 = 0. The general solution of (jHHJ) turns out to be 

V3/c 



x(t) 





(/ie- 


-kt _ 






c + 


m 




X 




/„'(' 






n 





cos(^t + 5)j"e^/o^' 



(4? 



where C3 = k 3 and Ji,C and 5 are arbitrary constants. By choosing n — 1, m = 
1, j(t) = and P(t) = k in equation (J3H|) one can get a special case of the Chazy 
equation XII (with N = 2 and parametric restrictions A = C3 and 5 = in Ref. [18]) 



x +4kxx + 3kx + 6k x x + k x + C3X = 0. 
The general solution of (}4*9*j) can be fixed easily from ()48|) in the form 



x(t) 



he~ kt + ei* cos(4^ + 5) 



C - ^(2/ie- fct - e3*(cos(*f*t + <*) + V3 sin(^* + 5)) 



(49) 



(50) 



3-4- Case (iv) q = constant, i = 1,2,3 

In this case the nonlocal transformation transforms the nonlinear ODE to the linear 
ODE (J35)) . The general solution of the latter becomes 



x(t) 



mit 



+ / 1 e m2i + / 2 e ra3f "ev/o^ 1 ' 



rn 



m\t' 



+ he m ^' + I 2 e 



m 3 t 



t \ n m 

e ™ 



(51) 
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where I\, I 2 and C are arbitrary constants and mf + Cim? + c 2 mj + c 3 = 0, i = 1,2, 3, 
and m 23 = |[— (mi + c\) ± \J—3ml — 2c\ui\ + c± — 4c2~]. 

Equation ()36)1 includes a class of important nonlinear systems. For n = 1, m — 
1, 7(f) = and (3{t) = k equation (J3T?)) becomes 

x +(ci + Akx)x + 3fcx 2 + (c 2 + 3A;ci + 6k 2 x)xx 

+ (ci + kx)k 2 x 3 + c 2 kx 2 + c 3 x = 0. (52) 

Equation (J52j) is a generalized version of the above special cases (JIHJ) and (}49|) of the 
Chazy equation XII, with additional terms. The general solution of (|52jl can be fixed 
easily from ()51|) in the form 



x{t) = [e mit + he m2t + J 2 e m3 ^ 



mi m 2 m3 



C + k[ - — + —e m2t + —e m ' il 



(53) 



To our knowledge, the above solution is a new one. For a detailed description of Chazy 
class of third order ordinary differential equations, their analytic properties and their 
general solutions one may refer to Ref. [18], besides the original work of Chazy [16]. 

In this section, we demonstrated the nonlocal connection that exists between linear 
and nonlinear ODEs of third order. In the following section, we extend the theory to 
nth order ODEs. 



4. nth order equations 



Let us consider a most general linear ODE of the form 

+ ... + C n _ 1 — + Cn)U(t) 



dt n + Cl dtl"-*) 



0, 



(54) 



where q's, i = 1,2, ...n, are arbitrary constants. The nonlocal transformation (fT^) 
connects (I54|) to the nonlinear ODE of the form 



+ C^ 1 ' + ■ ■ ■ + Cn-l^ + Cn )X 



Jn-1) 



(1) 



0, 



(55) 



where = (^ + (3(t)x m + 7(t)) n . For n = 2 and 3 equation (jHHjl coincides with 
(Unj) and (fBljj) respectively. Repeating the procedure described in Sees. 2 and 3 one can 
derive the general solution of the nth order nonlinear ODE, (j55)L in the form 



x(t) 



,£ Jo (a(f )-7(f))*' 



C + 



n jo 



/5(f)e^/o( a ('")-7(*"))^"^ t 



(56) 



Here again a(t) = -, where a is the known solution of equation (|54jl . The solution 
contains n + 1 integration constants, namely, Ji,/ 2 ,...J n and C. However, as was 
done in the second and third order ODEs, respectively, we can absorb one integration 
constant and rewrite the solution in such a way that it contains only n integration 
constants. Doing so we arrive at 

,a(t). 



x{t) = (— )»e n J o ,y ' 

J- n 



C + 



m 



n Jo V 4 



(57) 
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By appropriately fixing the constants and arbitrary functions in equation (J55|) one 
can deduce the nonlinear equation of interest. The general solution can also be fixed 
unambiguously from equation (f5Tj) . 



5. General Theory 

So far we focussed our attention only on the cases in which the nonlocal transformation 
that connectes the nonlinear ODEs with constant coefficient linear ODEs. However, one 
can also explore the hidden connection between nonlinear ODEs and variable coefficient 
linear ODEs. In the following we briefly discuss the underlying theory. 
Let us begin with a nth order linear ODE of the form 

/ d n S n ~^ d \ 

[dF + Cl{t) dK^) + • • • + Cn - l{t) Jt + c ^) u ^ = °. ( 58 ) 

where Cj's, % = 1, 2, . . . n, are functions of t, whose general solution be U = a(t), where 
a(t) is some known function. Now we consider a general nonlocal transformation of the 
form 

U{t) = g{t,x)eti f{t '' x)dt \ (59) 

where f(t,x) and g(t,x) are arbitrary functions of t and x, and substitute it into (|55|) 
so that the latter becomes a nonlinear ODE of the form 

Dt ] + c^Dt^ + ■■■ + c„-i(t)4 1} + Cn(t))g(t, x) = 0, (60) 

where = (4 + f(t,x)) n . Since we have assumed that the multiplicative function 
g(t,x) is a function of the variables t and x, the associated nonlinear ODE, ()60j) . takes 
a very complicated form. However, the task is to deduce the interesting cases from (|60J) 
which can be integrated explicitly. Again using the identity 

^ = M±i£ + / , (61 ) 
U g 

with U = a(t), equation (JBTj) can be brought to the form 

x = a t - f(t,x) 7—^, a=-. 62 

V / g x \t, x) g x (t,x) a 

Now, it is well known that only for certain specific forms of / and g, equation 

can be integrated [23]. At least for these cases, the general solution for the nonlinear 

ODE ()60|) can be obtained explicity. In this way one can classify a class of integrable 

nonlinear ODEs of any order through nonlocal transformations of the form (|59jl . 

We note that the equation (pjj) can also be obtained from equation by using a 

simple "gauge transformation" 

d n „ (n) (d „,. A n ... , d 



D 



dt n h 
and 

U(t)^g{x,t). 



(jt + f&^J' f{t,x) = -h{t,x) 
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Finally, we mention that in the case f(t,x) = the nonlocal transformation, (|59|) . 
becomes purely a local one and the resultant linearizable equations can be directly 
obtained from equation (jfiOj) . 

6. Conclusion 

In this paper, we have developed a novel way of identifying integrable nonlinear ODEs 
by connecting linear and nonlinear oscillator equations of any order through suitable 
nonlocal transformations. The proposed method is simple and straightforward. More 
importantly, our procedure not only identifies a class of integrable nonlinear ODEs of 
any order but also unambiguously gives their underlying solutions and thereby leads 
to the complete understanding of the dynamics of the given nonlinear system. As we 
have pointed out in the introduction, the modified Emden type equation (JTJ admits 
certain unusual nonlinear dynamical properties. The dynamical properties underlying 
other integrable nonlinear ODEs still remain to be explored in detail and are possibly 
worth investigating further. It is also of interest to investigate the existence of the above 
types of nonlocal connections in the case of partial differential equations as well, which 
we hope to pursue further. 
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